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ABSTRACT 

The lifetime of protoplanetary discs is intimately linked to the mechanism responsible 
for their dispersal. Since the formation of planets within a disc must operate within the time 
frame of disc dispersal it is crucial to establish what is the dominant process that disperses 
the gaseous component of discs around young stars. Planet formation itself as well as photo- 
evaporation by energetic radiation from the central young stellar object have been proposed 
as plausible dispersal mechanisms. There is however still no consensus as what the dominant 
process may be. In this paper we use the different metallicity dependance of X-ray photoe- 
vaporation and planet formation to discriminate between these two processes. We study the 
effects of metallicity, Z, on the dispersal timescale, fphot, in the context of a photoevaporation 
model, by means of detailed thermal calculations of a disc in hydrostatic equilibrium irradi- 
ated by EUV and X-ray radiation from the central source. Our models show fphot Z'^'^^ for 
a pure photoevaporation model. By means of analytical estimates we derive instead a much 
stronger negative power dependance on metallicity of the disc lifetime for a dispersal model 
based on planet formation. 

A census of disc fractions in lower metallicity regions should therefore be able to dis- 
tinguish between the two models. A recent study by Yasui et al. in low metallicity clusters 
of the extreme outer Galaxy ([O/H] — 0.7dex and dust to gas ratio of -0.001) provides pre- 
liminary observational evidence for shorter disc lifetimes at lower metallicities, in agreement 
with the predictions of a pure photoevaporation model. While we do not exclude that planet 
formation may indeed be the cause of some of the observed discs with inner holes, these ob- 
servational findings and the models and analysis presented in this work are consistent with 
X-ray photoevaporation as the dominant disc dispersal mechanism. 

We finally develop an analytical framework to study the effects of metallicity dependent 
photoevaporation on the formation of gas giants in the core accretion scenario. We show that 
accounting for this effect strengthens the conclusion that planet formation is favoured at higher 
metallicity. We find however that the metallicity dependance of photoevaporation only plays a 
secondary role in this scenario, with the strongest effect being the positive correlation between 
the rate of core formation and the density of solids in the disc. 

Key words: accretion, accretion discsxircumstellar matter- planetary systems:protoplanetary 
discs - stars ;pre-main sequence 



1 INTRODUCTION 

It is not yet established what is the dominant process that disperses 
the gaseous component of the discs around young stars. What is 
clear observationally is that evidence of gas accretion onto stars, to- 
gether with infrared/submm diagnostics associated with small dust 
grains entrained in the gas, disappear in young stars at an average 
age of a few Myr (Haisch, Lada & Lada 2001). This figure needs 
to be interpreted as an average figure since clearly individual stars 
may lose their discs on timescales that differ from this by at least 
a factor three (Armitage, Clarke & Palla 2003). It is equally clear, 
based on the relatively few systems observed in a state of transi- 



tion between disc possessing and discless status, that the timescale 
for disc dispersal (specifically the timescale over which regions of 
the disc are optically thin in the infrared) is much shorter than the 
overall disc lifetime (Skrutskie et al. 1990, Kenyon & Hartmann 
1995, Duvert et al. 2000). The fact that transition discs are relatively 
rare (constituting around 10% of the population of young stars with 
discs) is however an observational hindrance when it comes to es- 
tablishing the mechanism for disc clearing: transition discs appear 
to be a rather heterogeneous class of objects and it is hard to classify 
their diversity when the number of well studied objects in nearby 
star forming regions is still in single figures. An important feature 
of many transition discs, however, is that their spectral energy dis- 
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tributions are best fit by inner holes in the disc (these holes are not 
necessarily devoid of dust, but there is a large contrast in surface 
density between the outer disc and inner cavity) Q 

Although an obvious mechanism for the disappearance of dust 
diagnostics in discs is the simple coagulation of grains into entities 
that are large compared with the wavelength of observation (Dulle- 
mond & Dominik 2005), this scenario does not of itself explain 
why, in general, there is a correlation between the disappearance of 
dust and gas accretion diagnostics. Moreover, it is not obvious why 
dust coagulation should create a well defined inner hole structure 
rather than homogeneous depletion at all radii (although models 
combining grain growth with photophoresis - once the inner disc is 
optically thin to visible light in the radial direction - are promising 
in this regard; Krauss et al. 2007). 

The two leading disc dispersal mechanisms that satisfy these 
considerations are the formation of giant planets and photoevapo- 
ration. Both these mechanisms create an inner hole in the disc and 
also disrupt the accretion flow onto the star to some extenj^ In the 
case of photoevaporation, inner hole creation is succeeded by the 
rapid clearing of the outer disc thereafter. In the case of planet for- 
mation, clearing of the outer disc requires further planet formation 
at larger radii, i.e. the relative scarcity of transition discs argues that 
planet formation at one radius is rather rapidly followed by a wave 
of successive planet formation events in the outer disc ( Armitage & 
Hansen 1999). Some attempts have been made to classify whether 
individual transition discs are likely to be generated by planet for- 
mation or photoevaporation (Najita et al. 2007, Alexander & Ar- 
mitage 2007, Alexander 2008, Cieza et al. 2008, Kim et al. 2009) 
based on estimates of the accretion rate and disc mass. These analy- 
ses are however based on extreme ultraviolet (EUV) photoevapora- 
tion models (HoUenbach et al. 1994, Clarke et al. 2001, Alexander, 
Clarke & Pringle 2006a,b) which predict low accretion rates and 
disc masses in transition objects, in contrast to many of the ob- 
served systems. More recent models based on X-ray photoevapora- 
tion (Ercolano et al 2008, 2009, Owen et al. 2009) predict higher 
accretion rates and disc masses during transition and thus some- 
what muddy the observational distinction between the two classes 
of disc dispersal mechanism. Nevertheless it would appear clear 
that no single mechanism can explain all observed inner hole sys- 
tems: for example, photoevaporation is incompatible with the pre- 
dominantly high accretion rates in the sample of Kim et al. (2009) 
whereas a planet is perhaps an unlikely candidate for many of the 
objects contained in Cieza et al. (2008), which have very large holes 
and low accretion rates. 

At this point it is perhaps worth noting - despite the obvious 
interest of observed transition objects - that they do not necessar- 
ily hold the key to the process that disperses the majority of discs 
around young stars. This is because objects observed in transition 
can in principle be a mixture of stars that spend a significant time in 
this state (and must therefore be a minority of all young stars on sta- 
tistical grounds) and a general population of 'typical' young stars 
that must pass through such a transition quickly. (The argument 
that some fairly abrupt end to a disc's lifetime is required is well 

' This situation may not extend to discs around lower mass stars, as there 
are recent claims that discs in M stars may instead pass through an extended 
phase in which their dust content is homogeneously depleted at all radii 
(Currie et al. 2009). 

^ More or less completely in the case of photoevaporation (Clarke et al 
2001, Alexander, Clarke & Pringle 2006b) and to a variable extent in the 
case of planet formation, depending on the mass of the planet (Lubow et al. 
1999, Rice et al. 2003). 



illustrated by considering what would happen if a disc merely con- 
tinued to accrete onto the central star as a result of viscous evolution 
- the observed accretion rates, ages and disc masses of T Tauri stars 
would imply that integrated forward in time they would remain op- 
tically thick for > lOOMyr and then would linger for a comparable 
period as partially cleared systems). Therefore - regardless of what 
process explains the bulk of observed transition discs - we also need 
to decide what is the mechanism that terminates the lifetime of the 
majority of disc bearing young stars. 

In this paper we therefore take a different approach to dis- 
criminating between planet formation and photoevaporation as the 
prime disc dispersal mechanism, by instead considering the metal- 
licity dependence of these two processes. Although most studies of 
circumstellar discs have currently been undertaken in regions with 
a rather low range of metallicities, the advent of the Herschel Space 
Observatory and later the James Webb Space Telescope raises the 
prospect of being able to conduct disc censuses in more distant re- 
gions whose metallicity differs considerably from that of star form- 
ing regions in the solar neighbourhood. Such studies are already 
possible in the Extreme Outer Galaxy (EOG) with current instru- 
mentation as demonstrated by the recent work of Yasui et al. (2009) 
which we will discuss in more detail later. Our focus here will not 
be on the nature of the (in any case rare) transition objects, but will 
instead centre on how the mean disc lifetime is expected to vary 
with metallicity in these two scenarios. To this end we combine 
simple models for disc viscous evolution with (metallicity depen- 
dent) photoevaporation and also with semi-analytic prescriptions 
for planet formation (Section 2). We will show that in the photoe- 
vaporation model the disc lifetime is a mildly increasing function of 
metallicity, resulting from the rather higher photoevaporation rates 
in the case of low metallicity gas for which opacities are lower and 
line cooling less efficient. On the other hand, the disc lifetime is 
a strongly decreasing function of metallicity in the case that disc 
clearing is initiated by planet formation. This is simply because 
at higher metallicity the assumed higher surface density of plan- 
etesimals in the disc encourages the more rapid formation of a gas 
giant planet, specifically because it accelerates the formation of a 
rocky core of the critical mass required to initiate the accretion of 
a gaseous envelope. This is borne out, for example, by the hydro- 
dynamical models of gas giant formation of Pollack et al. (1996) 
and Hubickyj et al. (2005lj. This positive link between metallic- 
ity and planet formation has been noted by a number of authors 
(e.g. Ida & Lin 2004b) as a possible explanation of the observed 
increase in the frequency of giant planets as a function of metallic- 
ity (Santos, Israelian & Mayor 2000, 2001, 2004; Gonzalez et al. 
2001 ; Sadakane et al. 2002; Heiter & Luck 2003; Laws et al. 2003; 
Fischer & Valenti 2005; Santos et al. 2002 and references therein) 
but has not been previously linked to the issue of disc lifetimes. In 
Section 3, we revisit the connection between the frequency of giant 
planets and metallicity in the case of a hybrid model where planet 
formation takes place in the context of a nebula that is subject to 
metallicity dependent photoevaporation of the disc gas. Since pho- 
toevaporation is more effective at low metallicities and thus reduces 
the disc lifetime, this strengthens the conclusion that planet for- 
mation is favoured at high metallicity. Nevertheless, the additional 



^ Note that this conclusion is based on an assumed linear relationship 
between the density of planetesimals and the metallicity; this conclusion 
would be only strengthened if one takes into account the recent suggestion 
of Anders et al 2009 that the efficiency of planetesimal formation by the 
streaming instability should increase steeply at higher metallicity. 
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Figure 1. The metallicity dependance of mass loss rates (left panel) and dispersal timescales (right panel) of protoplanetary discs photoevaporated by X-ray 
and EUV from the central star. Right Model data points (asterisks) obtained from our photoionisation modeling; the dashed line is a least square fit to these 
data. Left: The asterisks show the photoevapoaration timescales obtained substituting the points on the right panel into equation AlO and normalising to 2 Myr 
for solar metalHcities. The dashed line shows the relation given in equation 4 for p = 1. 



ingredient of metallicity dependent photoevaporation is found to 
be a second order effect in determining the positive correlation be- 
tween planet frequency and metallicity. We present our conclusions 
in Section 4. 



2 THE DEPENDENCE OF DISC LIFETIME ON 
METALLICITY 

2.1 Lifetime against photoevaporation 

We consider the case where the gas in the disc (initially with mass 
Mgo) undergoes viscous accretion onto the central star. As a re- 
sult of this secular evolution, the disc radius grows while the disc 
mass and the accretion rate onto the star decline with time. We here 
parameterise this process according to the similarity solutions that 
pertain in the case that the kinematic viscosity in the disc is a sim- 
ple power law function of radius (y oc r''; Lynden-Bell & Pringle 
1974, Hartmann et al 1998). This is not likely to be a good descrip- 
tion in detail since even in a simple a disc model the viscosity is a 
function of both surface density and radius, depending on the dom- 
inant local opacity source; yet further complications are introduced 
when more realistic prescriptions for angular momentum transport 
- involving the action of the magnetorotational instability or of self- 
gravitating torques are included (e.g. Armitage et al. 2001, Clarke 
2009, Rice & Armitage 2009, Zhu et al 2009, Cossins et al. 2009). 
The prescription we use here however has the advantage of provid- 
ing a simple analytic form such that the value of p is linked to the 
power law index of the surface density profile (i.e. over much of the 
radial range of the similarity solution, the disc's surface density is 
a power law E oc r"''). Our choice of p (which fixes the power law 
exponents for the evolution of disc quantities) is therefore obser- 
vationally motivated by fits to the infrared spectral/submm energy 
distributions of young stars, which suggest p ~ \ (Andrews et al 
2009, Kitamura et al 2002, Isella et al 2009) 

The Appendix contains a heuristic derivation of the asymp- 
totic form of the similarity solutions; according to eguation lASI the 
accretion rate declines as; 

M oc Myj^h'-'^'-"^^^-^-^"^ (1) 



where rja is the initial disc scaling radius. Such a power law 
decline obviously does not enable one to define a disc lifetime, 
since the timescale on which disc properties change by order unity 
is always of order the present age. We therefore need to intro- 
duce a further condition for disc dispersal. In the case that we in- 
voke (metallicity dependent) photoevaporation, at a constant rate 
M = Mm/(Z), then the disc disperses rapidly at the point that the ac- 
cretion rate through the disc falls to Miy(Zl]. Thus the disc lifetime 
against photoevaporation can be written; 

t,„„ oc Mv.(Z)<*-2">/<-'-2'"M<:-^'"«^-^">r<^-''«'-^''> (2) 

At this point it becomes necessary to specify the source of the 
photoevaporative flow. Until very recently it has been widely as- 
sumed that such a flow is mainly driven by the extreme ultraviolet 
(EUV) photons from the central star However we have recently 
shown that soft (0. 1-lkeV) X-ray irradiation can drive powerful 
winds with typical mass loss rates of 10"'-10"* Mg/yr (Ercolano 
et al., 2008b, 2009 ECD09, Owen et al., 2009). In the hydrostatic 
equilibrium model of ECD09, X-rays from the central pre-main- 
sequence star ionise and heat the gas in the disc atmosphere; the 
photoevaporation rate is estimated by assuming that gas that be- 
comes hotter than the local escape temperature becomes unbound 
and escapes at the sound speed. In reality, pressure gradients within 
the disc are such that a photoevaporative flow can in fact be initi- 
ated subsonically in deeper layers where the temperature is less 
than the escape speed. This produces even larger winds (Owen et 
al. 2009) which exceed photoevaporation rates from EUV radiation 
by two orders of magnitude. In fact, for strong X-ray emitters (i.e. 
with Lx ~ 10^° erg s"'), these large flow rates are similar to the me- 
dian accretion rates measured in T Tauri stars (Natta et al. 2006). 



Here we are ignoring the fact that, as shown by Owen et al. (2009) the 
rapid disc dispersal begins once the accretion rate has fallen to roughly ten 
times below the photoevaporation rate, causing the system to experience 
a period of photoevaporation starved accretion (Drake et al. 2009) prior 
to the fast dispersal phase. These consideration however do not affect the 
proportionality relations and therefore the conclusions of this work remain 
unchanged. 
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When one also takes account of the dispersion in X-ray luminosi- 
ties among young stars (Preibisch et al. 2005, Albacete-Colombo 
et al. 2007) and the associated range of expected photoevaporation 
rates, it then becomes plausible that - in the absence of competing 
disc dispersal mechanisms - the lifetime of gas discs in T Tauri stars 
is terminated by X-ray powered photoevaporation. 

With this in mind, we therefore use the metallicity dependence 
of the photoevaporation rate derived from our X-ray photoevapo- 
ration models. These employ the modeling strategy described by 
Ercolano et al. (2008b) and ECD09 to obtain the temperature struc- 
ture and photoevaporation rates for gaseous discs with a range of 
metallicities, spanning from 0.01 solar to twice solar. Briefly, this 
involves the coupling of the 3D photoionisation and dust radia- 
tive transfer code MOCASSIN (Ercolano et al. 2003, 2005, 2006, 
2008a) with an iterative solution of the hydrostatic equilibrium 
structure of the irradiated disc. Our initial guess at the density dis- 
tribution is the disc model of d'Alessio et al (2001) for a 4000 K 
0.7 Mo pre-main sequence star surrounded by an optically thick 
flared circumstellar disc with an outer radius of 500 AU and a total 
mass of 0.027 M©. This structure is irradiated by an unscreened 
EUV-l-X-ray spectrum (see description in ECD09). We refer to 
ECD09 for further details including the dust model and the irra- 
diating spectrum. Our solar abundance set is that of Asplund et 
al. (2005) depleted according to Savage & Sembach (1996). The 
dust to gas ratio for the solar abundance case was set to 6.5x10"'' 
(see discussion in ECD09). For models of different metallicity, the 
dust to gas ratio and the abundances of metals were multiplied by 
(Z/Zq). As noted above, we estimate the photoevaporation rates by 
assuming that a sonic flow is launched from regions where the tem- 
perature exceeds the local escape temperature. Our preliminary hy- 
drodynamic calculations (Owen et al. 2009) suggest that this may 
under-estimate the total mass flow rates, although we do not expect 
this to have a large effect on the relative scaling with metallicity. 

The left panel of Figure [T] shows the dependance of the to- 
tal photoevaporation rates on metallicity, Mw(Z), which can be 
approximated by a power-law of index -0.77. The model data 
points shown in the left panel of this figure were obtained using 
the approach described above, which involved rerunning model 
FS0H2Lxl from the set described by ECD09 for the appropriate 
metallicities and dust-to-gas ratios. The corresponding photoevap- 
oaration timescales, shown in the right panel, were then obtained 



by substituting these numerical values of Mw(Z) in equation A 10. 
The resulting timescales were then normalised such that a lifetime 
of 2 Myr is obtained at solar metallicities. 

The increasing X-ray photoevaporation rates at lower metal- 
licity can be readily understood by looking at the gas temperature 
distributions in the disc for the different metallicity cases that are 
shown in Figure |2] The reduced extinction in the low metallicity 
cases allows high density gas at larger columns to be ionised and 
heated to temperatures sufficiently high for the gas to be entrained 
into a photoevaporative flow. Other metallicity-dependent efl'ects 
are also present and work in the same direction (i.e. lower metal- 
licity = higher gas temperatures), but they play a secondary role. 
These efl'ects are listed here for completeness: a low metallicity im- 
plies a reduced cooling by fine structure lines of ions and neutrals 
(such as [O i] and [C ii], e.g. Ercolano et al. 2008b), which leads 
to higher temperatures in the X-ray heated disc atmosphere, as par- 
ticularly evident in the right panel of Figure |2] Also a lower dust 
to gas ratio means that the gas experiences a reduced competition 
from the dust grains for the absorption of the energetic photons and 
that the dust-gas collisional cooling term, which may provide a sig- 
nificant contribution at the base of the photoevaporative envelope, 
is reduced at lower metallicities. We emphasise however that by 
far the dominant mechanism here is the reduced opacity which in- 
creases the penetration column of the ionising radiation. It is worth 
noting at this point that such a large effect is not expected to occur 
for gas ionised by EUV only. X-ray photons are mainly absorbed 
by the inner shells of the more abundant heavy elements in the gas 
and dust (e.g. oxygen, carbon etc.) while for EUV photons by far 
the largest source of opacity is hydrogen, which is not affected by 
changes in metallicity (although a reduction in the dust-to-gas ratio 
would still reduce the total EUV opacity). We also note that a pho- 
toevaporation process based on far-ultraviolet (FUV) irradiation, 
would also follow a different dependance on metallicity (Gorti & 
Hollenbach 2009). Therefore the predicted metallicity dependance 
shown here pertains only to X-ray photoevaporated discs. 

When this metallicity dependent wind mass loss rate is com- 
bined with equation lAlOl one obtains the result that 

f nr 70.77(4-2p)/(5-2p).,(4-2p)/(5-2p) (2-p)/(5-2p) 

'pilot ''igQ r^Q (J) 

Since we do not expect the distribution of initial disc gas prop- 
erties (Mgo and r^g) to depend on metallicity, we thus deduce that 
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the distribution of disc lifetimes in regions of different metallicity 
should scale with a Z dependent factor, i.e. 



t nr 70.77(4-2p)/(5-2p) 



(4) 



This Z dependence is rather weak, with a power law exponent 
of 0.52 (for p = I) and 0.38 (for p = 1.5). The right hand panel of 
Figure [T] is an illustrative example of the dependence on metallic- 
ity of the mean disc lifetime of a population of discs (with p = I) 
that is normalised such that the mean disc lifetime for discs of so- 
lar metallicity is 2Myr. We emphasise that this normalisation is set 
by our assumptions about the disc's secular evolution timescale, 
as controlled by its equivalent viscous 'a' parameter (Shakura & 
Sunyaev 1973). We can readily adjust this quantity within a plausi- 
ble range in order to reproduce an observationally reasonable mean 
disc lifetime. An increase in metallicity from solar to twice solar 
causes the disc lifetime to increase from ~2 to ~3.1 Myr, and a 
decrease from solar to -0.7 dex solar, as in Cloud 2 in the EOG 
observations of Yasui et al. 2009, produces a decrease in the disc 
lifetime from ~2 to -0.7 Myr, which is roughly consistent with the 
observed disc fractions. 

2.2 Disc lifetime against gas giant planet formation 

We now assume that disc dispersal is instead initiated by the for- 
mation of the first gas giant planet, as is a common interpretation 
of observed transition discs (Rice et al. 2003, Setiawan et al. 2007). 
According to the core accretion model, the formation of a gas giant 
involves first the accumulation of a solid core of critical mass (of 
order 10-20 earth masses) followed by the hydrodynamic accretion 
of gas from the surrounding disc (Pollack et al. 1996). This latter 
occurs on the Kelvin Helmholtz timescale of the planet, which is 
~ a Myr for a 20 earth mass planet and decreases strongly with 
increasing planet mass. For the purpose of our analytic estimates 
(following Ida & Lin 2004a) we therefore assume that the limiting 
timescale in the creation of a gas giant planet (i.e. a roughly Jupiter 
mass planet which can clear a gap in the disc) is the time required 
to accumulate a core containing a critical mass of solids. In fact hy- 
drodynamic calculations of planet formation paint a more complex 
picture in which the majority of the formation process is instead 
spent in a phase ('Phase 2';Hubickyj et al. 2005) where the planet 
accretes a mixture of gas and planetesimals and where, depending 
on the properties of the background nebula, gas giants can form, 
albeit more slowly, with a core mass somewhat below the critical 
value mentioned above (we discuss below how such considerations 
alfect our conclusions). 

Proceeding for now with our simple analytic estimate we eval- 
uate the timescale needed to achieve a fixed critical core mass (to- 
gether with its scaling with metallicity), by using the analytic ex- 
pression for the growth of solid core mass (M^) given in Ida & Lin 
(2004a, IL04): 



dt 



(5) 



where is the surface density of rocky planetesimals, and Sg 
is the instantaneous local value of the gas surface density. We fol- 
low IL04 by assuming that the planetesimal distribution is decou- 
pled from the evolutionary processes that control the evolution of 
Eg. We thus assume that a fixed fraction of the initial solid content 
of the disc forms a planetesimal disc whose radial profile follows 
the initial radial profile of the gas and that is then constant in 
time (we thus neglect the depletion in due to the accretion of 
planetesimals during core growth; this is acceptable to first order 



provided that the core is not close to achieving its isolation mass, 
an issue to which we return below.) 

We however (see Appendix) differ from IL04 in that we use 
the viscous similarity solutions to determine the time dependence 
of Sg rather than adopting an ad hoc exponential reduction in the 
gas surface density. Solving equation IA13l and requiring that the 
core mass attains a fixed value for runaway gas accretion (see Ap- 
pendix) then implies a formation timescale: 



^form ^ Z 



S(2-l>)/{S-3p)j^-n2-p)H5-3p)J2-p}(9-5p)/(S-3p)j2-p)ap+3)/(5-3p), 



(6) 



Since this timescale is an increasing function of radius, it fol- 
lows that planet formation is first favoured at small radius. However 
this effect is counteracted by two effects that favour core growth 
at larger radius. Firstly, the sticking efficiency of planetesimals in- 
creases once the radius at which ice sublimes (a,„) is exceeded 
(Hayashi 1981, Pollack et al 1994). Secondly, it is necessary that a 
critical core mass can form without consuming all the planetesimals 
in its 'feeding zone' . This latter is often assumed to be an annular 
region around the core with width related to the Hill radius of the 
core (i.e. of width that scales linearly with orbital radius). Thus the 
radius {a,g) at which a core of critical mass just consumes all the 
planetesimals in its feeding zone obeys a scaling of the fornjf] 



Erffl, = constant 



«„<xZ-"<2-"V„oM-'«^-'" 



(7) 



(8) 



Therefore if > i^tg (i-e- if the minimum radius for planet 
formation is set by ice sublimation rather than feeding zone consid- 
erations) we have: 

tf...,..,„o.iz-^Mj,^^^f-'''''-'^' (9) 
Alternatively, if a,g > aice we substitute[8]into[6land obtain 



'■form,jj 



:(z(-13-2p)^^-n^(4-2,„(6+„)y/(5-3/') 



(10) 



We thus see that these planet formation timescales are highly 
Z dependent, in the sense that planet formation is faster in discs 
with higher metallicity. For p = I and p = 1.5 equation |9] implies 
that the planet formation timescale scales either with Z"^^ or Z^^, 
while[TO]implies scaling either with Z"'-'' or Z"'- 

This extreme level of Z dependence is almost certainly an ar- 
tifact of the fact that we have assumed that the timescale for giant 
planet formation is completely controlled by the timescale for solid 
core growth. It is also worth mentioning that this very steep scaling 
stems from the fact that the prescription for core growth (Ida & Lin 
2004a; equation 5 above) involves a dependence on gas density; 
physically, this accounts for the fact that gas drag reduces the ve- 
locity dispersion of the planetesimals and hence increases the core's 
accretion cross section due to gravitational focusing (Kokubo & Ida 



^ If, as discussed in Ida & Lin (2004), one relaxes this assumption so that 
the core can feed from the entire stock of planetesimals interior to its orbit 
then one also obtains a scaling of the same form. 

* We also note that the planet formation timescale is highly dependent on 
the initial disc mass, implying that planet formation should only be possi- 
ble in systems belonging to the upper centiles of the disc mass distribution. 
When the disc mass threshold implied by the Ida & Lin models is combined 
with observations of the disc mass distribution by Andrews & Williams 
(2005), the predicted planet frequency is consistent with those found in ra- 
dial velocity surveys at solar metallicity: see discussion in Wyatt, Clarke & 
Greaves (2007). 
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2002). We assume that the disc gas surface density declines as a 
power law due to viscous evolution and thus the core grows as a 
low power of time. Thus the time required to achieve a given core 
mass is a strong function of the normalisation of the solid mass. 
This is not necessarily incorrect and is in this respect probably more 
realistic than models of core accretion which adopt a temporally 
constant background gas surface density (to date such models ei- 
ther keep this quantity constant or introduce a linear decline to zero 
over 1-3 Myr; Lissauer et al. 2009). 

We now turn to the results of hydrodynamic models of planet 
formation and see to what degree the limited range of available 
models in fact support our estimates above. These models demon- 
strate two effects: increasing the surface density of solids decreases 
the timescale for planet formation, in qualitative agreement with the 
arguments above. Another factor in the opposite direction is that at 
higher metallicity the opacity of the accreting gaseous envelope is 
increased; this slows the contraction of the envelope and thus slows 
down gas accretion. (Note that our simple model above only con- 
siders the time needed to acquire a critical solid core mass and does 
not consider the subsequent time spent in acquiring a similar mass 
in gas which is considerable in the hydrodynamic calculations; Hu- 
bickyj et al. 2005). Nevertheless, it is found that the acceleration of 
planet formation due to enhanced planetesimal surface density is a 
much stronger effect than the retarding effect of increased opacity. 
For example, Hubickyj et al. found that the time to form a gas gi- 
ant increases by around a factor three if the dust opacity is boosted 
by a factor fifty (see also Ayliffe & Bate 2009); on the other hand 
a modest increase of the planetesimals surface density by a fac- 
tor 1.6 reduces the formation timescale by a factor six. Thus the 
simulations support our qualitative conclusion that the net effect 
of increasing the metallicity is to strongly decrease the formation 
timescale of giant planets 

More detailed comparison is not really warranted given the 
fact that the choice of simulations performed was motivated by 
different scientific questions than the one we ask here. We are in- 
terested in the time required to form a planet at given metallicity 
anywhere in the disc. Since we have imposed a hard requirement 
that the growing protoplanet should achieve a fixed core mass, we 
exclude planet formation in inner regions of the disc where the iso- 
lation mass is less than this critical value; in our prescription we 
first form a planet at a larger radius where the formation timescale 
is longer. Evidently this slows down planet formation at low metal- 
licity. The hydrodynamical calculations are instead motivated by 
creating Jupiter at a particular radial location. In low metallicity 
runs planet formation takes longer because the forming protoplanet 
has to 'make do' with a lower core mass and this then greatly in- 
creases the timescale for accreting a comparable mass of gas. In a 
real disc both these effects would come into play and it is not clear 
whether, under these circumstances, the first planet to form is at 
larger radius (where the core mass is higher but the core accumula- 
tion timescale is longer) or at smaller radius where the bottleneck is 
instead the slow accumulation of gas onto a smaller mass rock core. 
What is important to our discussion is that both effects are in the 
same direction and imply a strong (negative) dependence of planet 
formation timescale on planetesimal surface density (and thus im- 
plicitly metallicity). 

If the majority of protoplanetary discs dispersed due to planet 
formation then a disc census in a lower metallicity region, such as 
Cloud 2 in the EOG, would yield much longer disc lifetimes than 
those derived in the solar neighbourhood. This is in contradiction 
with the recent observations of Yasui et al. (2009) who find disc 



lifetimes shorter than 1 Myr, compared to the few Myr found in the 
solar neighbourhood (e.g. Haisch et al. 2001). 



3 THE METALLICITY DEPENDENCE OF PLANET 
FORMATION IN DISCS SUBJECT TO 
PHOTOEVAPORATION 

As the list of exoplanets discovered in the solar neighbourhood con- 
tinues to grow, statistical significance is lent to the observation that 
solar type stars hosting giant planets are on average more metal-rich 
than field stars (Santos, Israelian & Mayor 2000, 2001, 2004; Gon- 
zalez et al. 2001; Sadakane et al. 2002; Heiter & Luck 2003; Laws 
et al. 2003; Fischer & Valenti 2005; Santos et al. 2002 and refer- 
ences therein). While the origin of this metallicity excess has been 
the object of much debate in the literature, with the two main sce- 
narios being "primordial" or "external" enrichment, there is a grow- 
ing consensus that metal rich gaseous discs are a more favourable 
environment for the formation of the rocky cores required by core 
accretion models (e.g. Ida & Lin, 2004a). 

The phenomenological simulations of gas giant planet for- 
mation of IL04 (whose assumptions imply, as we showed above, 
a strong reduction in planet formation timescale with metallicity) 
demonstrated a statistical preference for planet formation and sur- 
vival in discs with a high ratio of solids to gas. If one ignores the 
possibility of stellar surface contamination then the solid to gas ra- 
tio simply scales with the metallicity of the central star. Such sim- 
ulations therefore reproduce the observed positive correlation be- 
tween the incidence of planets and stellar metallicity. Indeed, to ze- 
roth order, the result results of IL04 can be understood by adopting 
their prescription for core growth (equation lAl lb and deriving the 
region of parameter space of initial conditions, at given Z, that lead 
to the creation of a critical mass core. Wyatt, Clarke & Greaves 
(2007) showed that this is roughly equivalent to thresholding the 
disc mass distribution at a fixed mass in solids; they showed that, 
using the observed disc mass distribution (Andrews & Williams 
2005) one can reproduce the observed dependence of planet fre- 
quency on metallicity under this simple assumption. 

As noted above, the models of IL04 omitted viscous evolution 
of the disc gas and instead depleted the gas on an adjustable global 
e-folding timescale. Moreover, they included no hard cut-off in disc 
gas lifetimes as would result from photoevaporation. We here there- 
fore revisit this issue under the assumption that the lifetime of the 
disc gas is terminated by (metallicity dependent) photoevaporation. 
In the context of our analysis in Section 2, we can readily establish 
which sets of disc initial parameters should lead to giant planet for- 
mation by looking at the subset of models for which ^/-onn,,,,,, < tpimt 
(equations[3l[9l[T0t. We then obtain the criteria: 

M,o > ■ ,2<2-rt'/(«-4">Z-<5-2''>/<»-^'')M^(z)2(5-3/')/5(9-4rt (j j) 

if a,„ > a,g and 

M,o>Kr{r-Jt'''''-'''''*''"'Z-''^^^^^^^^^ 

if a,„ < a,^. Ki and K2 are constants. 

We thus obtain the expected result that a higher threshold gas 
mass is required at low metallicities; therefore if we make the rea- 
sonable assumption that the disc gas mass distribution should be 
independent of metallicity one recovers the qualitative result that 
the incidence of giant planets should increase with metallicity. In 
order to quantify this, one of course has to make specific assump- 
tions about the form of the distributions of initial disc masses and 
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radii (Wyatt et al. 2007). Here, however, our interest lies elsewhere: 
we just want to discover whether the particular aspect of metallic- 
ity dependent photoevaporation rates is an important factor in driv- 
ing this relationship. To this end we compare the powers of Z in 
equationsll lland[T2lthat would result from the case that the photo- 
evaporation rate was independent of Z (i.e. Mw{Z) = constant) with 
the power of Z in the case that we employ our X-ray derived value 
M„iZ) oc Z""-^'. For p = 1 we find that this changes the power of 
Z in equation [TTI from -0.6 to -0.88 and in equation [T2l from -0.8 
to -0.83. This therefore implies that the metallicity dependence of 
the photoevaporation rate has a rather minor part to play in explain- 
ing the observed positive correlation between planet frequency and 
stellar metallicity. This can be readily understood inasmuch as we 
have already seen (equation I A lOt that the disc lifetime is not ex- 
pected to be a strong function of metallicity in this case and thus 
this does not in itself have a large impact on the planet producing 
capacity of a disc. 



4 CONCLUSIONS 

We have shown that the study of disc lifetimes in regions of dif- 
ferent metallicities can be used as a powerful discriminant between 
the two currently leading models of disc dispersal - photoevapora- 
tion and planet formation. By means of detailed thermal and pho- 
toionisation calculations we have determined that a disc's lifetime 
against photoevaporation is a shallow positive power of metallicity, 
fphot °^ Z"-^^. This metallicity dependance is specific to a photoevap- 
oration mechanism driven mainly by X-ray radiation as it relies on 
the significant reduction of the gas opacities with metallicity. We 
have also shown that, on the contrary, a disc dispersal mechanism 
based on planet formation yields disc lifetimes that are a strong 
negative power of Z. Therefore, a census of disc fractions in regions 
of lower metallicities compared to the solar neighbourhood will be 
crucial to determine which is the dominant mechanism responsible 
for the rapid demise of protoplanetary discs. 

Recent observations of embedded clusters in Cloud 2 of the 
Extreme Outer Galaxy (EOG) presented by Yasui et al. (2009, 
YS09) indeed find shorter disc lifetimes for this metal poor en- 
vironment ([O/H] — 0.7 dex), compared to the solar neighbour- 
hood. These results are in agreement with the predictions of an 
X-ray-l-EUV photoevaporation model, and argue against planet for- 
mation as the dominant dispersal mechanism. 

YS09 quote a mass detection limit of 0. 1 Mq for Cloud2 in the 
EOG using the 8.2m Subaru telescope, similar to values obtained 
for embedded clusters in the solar neighborhood with smaller (2- 
4m class) telescopes, and thus argue that a comparison of the disc 
fractions in embedded clusters in the EOG with those of clusters 
in the solar neighbourhood is justified. Further observations aimed 
at determining disc fractions in low metallicity regions of various 
ages would be very useful to confirm these important results. 

We finally show that, in the context of giant planet formation 
in the core accretion scenario, the effect of metallicity dependent 
photoevaporation is to strengthen the conclusion that planet forma- 
tion is favoured in high metallicity environments since the lifetime 
of the disc against photoevaporation (fphot) is a positive function of 
Z. This effect, however, only plays a secondary role: the main pro- 
cess that favours planet formation at high metallicity is simply the 
faster core growth in the case of a high surface density of solids in 
the disc. 
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APPENDIX A: ANALYTICAL DERIVATIONS 

In order to evaluate the metallicity dependence of giant planet for- 
mation, we need to establish how the background nebula evolves as 
a function of the initial conditions, i.e. the initial total disc mass in 
gas, Mgo, and the initial disc outer radius, rjQ. We assume that the 
evolution of the disc gas is governed by viscous evolution provided 
that the accretion rate through the disc exceeds the photoevapora- 
tion rate, Myf,{Z), and that at this point the disc gas is dispersed, 
thus ruling out subsequent gas giant planet formation. We start by 
providing a heuristic derivation of the asymptotic evolution of a vis- 
cous disc (see Lynden-Bell & Pringle 1974, Hartmann et al 1998) 
in which the kinematic viscosity is a power law function of radius, 
i.e. 

veer" (Al) 

We first note that such a disc evolves through a sequence of 
discs with increasing radius, rj, and decreasing gas mass, Mg, such 
that the total disc angular momentum is to first order conserved 
(since the fraction of the total disc angular momentum that is ad- 
vected onto the star is small). Thus, since most of the disc's angular 
momentum resides at large radius, we have: 



1/2 



(A2) 



Furthermore, the disc evolves on the viscous timescale at r^; 
since this is long compared with the viscous timescale at radii << 
rj, it follows that the disc at radii << is in an approximately 
steady state, i.e. the accretion rate, M, is independent of r. In this 
case, accretion disc theory (e.g. Pringle 1981) relates this (at radii 
well away from the disc's inner edge) to the kinematic viscosity 
and surface density, Sj(r), via 



(A3) 



It thus follows that the asymptotic form of the viscously evolv- 
ing disc over much of its radial extent is in this case given by the 
power law Eg oc r If we normalise this power law by the instan- 
taneous values of the disc mass and radius we then obtain 



^ "'-g'd 



(A4) 



Combining this with equation IA2I we can eliminate rj and 
obtain: 



(A5) 

Now from equations I A 1 1 1 A3 l and l A5 1 we can write 

from which we obtain the asymptotic power law scalings: 

Mg cx Mgor^/^r'/**-^'" (A7) 



M oc Mgori/'f<-^+2">/<^-^'" 
and 

Eg cc Mgor]^-r-"/<-^+2'"/(^-2"> 



(A8) 



(A9) 



If we now define the disc lifetime against dispersal by photo- 
evaporation as the disc age such that M = Mw{Z), we have: 



(AlO) 



We now consider planet formation in the context of such a 
viscously evolving gas disc and follow IL04 by assuming that the 
growth of the rock core is governed by an equation of the form 
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^^^,(r)r-"'M^^'l.f (All) 

Erf is the surface density of rocky planetesimals. We follow 
IL04 by assuming that this distribution is decoupled from the evo- 
lutionary processes that control the evolution of Eg. We thus assume 
that a fixed fraction of the initial solid content of the disc forms a 
disc whose radial profile follows the initial radial profile of the gas 
and that is then constant in time (we thus neglect the depletion 
in Srf due to the accretion of planetesimals during core growth; this 
is acceptable to first order provided that the core is not close to 
achieving its isolation mass, an issue that we consider further be- 
low). Thus 

Erf oc ZMjor^<'""V-'' (A12) 
substituting from lA9l and lA12h nto lAl ll we obtain: 

^ cc ZM;/V,^''-^""V-<7"+3'/5M2/3r<2/5)<5-2rt/(4-2„) (^13) 

We conclude from lA13] that the timescale for the formation of 
a rocky core of critical mass for the accretion of a gaseous envelope 
then scales with radius and initial disc parameters according to: 

tf„m OC Z"^'^<^"''''<^''"^^'''M"J''"''''*^"^''V*q"'''*''"^''''*^'^"'''V'-"'''*''''*^''<'"^'''(A14) 

Equation IA14I demonstrates that the formation timescale in- 
creases with increasing radius, as expected given the lower surface 
density and larger orbital timescale at larger radius. Thus planet for- 
mation occurs first at the minimum radius that is allowed according 
to two further criteria. First of all, it is necessary that this radius is 
at least as large as the ice sublimation radius (a,„) since the stick- 
ing efficiency of planetesimals is much lower in regions devoid of 
solid ice. Secondly, it is necessary that a critical core mass can form 
without consuming all the planetesimals in its 'feeding zone'. This 
latter is a region whose width is related to the Hill radius of the core, 
and thus is a linear function of orbital radius. Thus the radius (a,g) 
at which a core of critical mass just consumes all the planetesimals 
in its feeding zone obeys a scaling of the form 

Erffl^j = constant (A 15) 

i.e. 

fl„cxZ-'/P-"VrfoM;„'«^-''> (A16) 

Therefore if aj„ > a,g (i.e. if the minimum radius for planet 
formation is set by ice sublimation rather than feeding zone consid- 
erations) we have: 

fform,,,,,, 'X [Z M^^r^g "j (A 17) 

Alternatively, if a,g > a,„ we substitute IA16I into IA14I and 
obtain 

?,„„„,„„, cc (z<-'^-^'"m;o'V^:-^"><^^">)'"'-''" (A18) 

Finally, the condition that at least one gas giant planet is able 
to form prior to disc photoevaporation is given by the condition 

^/o™,„,„ ^ tphotf I.e. 

Mgo > Ki ■ ;-2P-")'/<''-^'')z-(5-2p)/(9-4rt^^(2)2<^^-3">/5(9-4p) (A19) 
if ajcc > a,g and 

if ai„ < a,g. Ki and K2 are constants. 
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